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We investigate weak nuclear decays with extremely small kinetic energy release (Q value) 
and thus extremely good sensitivity to the absolute neutrino mass scale. In particular, we 
£N| i consider decays into excited daughter states, and we show that partial ionization of the parent 

| atom can help to tune Q values to 1 keV. We discuss several candidate isotopes undergoing 

Qh' P^, bound state (3, or electron capture decay, and come to the conclusion that a neutrino 

' mass measurement using low-Q decays might only be feasible if no ionization is required, 

and if future improvements in isotope production technology, nuclear mass spectroscopy, and 
^vq | atomic structure calculations are possible. Experiments using ions, however, are extremely 

challenging due to the large number of ions that must be stored. New precision data on 
nuclear excitation levels could help to identify further isotopes with low-Q decay modes and 
O ■ possibly less challenging requirements. 

PACS numbers: 14.60.Pq, 23.40.-s, 23.40.Bw 
> ■ 1. INTRODUCTION 

o\ 

One of the big unknowns in astroparticle physics today is the absolute neutrino mass scale 
m v . While indirect probes such as cosmology [l[ and neutrinoless double /3 decay 0] achieve sub- 
eV sensitivity to m v , it is desirable to complement these measurements with model-independent 
direct bounds. The most advanced efforts in this direction have been the kinematical studies of 
the j3 spectrum in tritium decay by the Mainz U and Troitsk 0] collaborations, yielding the limit 
m u < 2 eV. In the near future, the sensitivity will be improved to m u < 0.2 eV by the KATRIN 
experiment [fjj]. However, Mainz, Troitsk, and KATRIN are limited by the accuracy to which 
the spectrum of decay electrons can be measured few eV below the kinematical endpoint, where 
the impact of m u > is largest. Since the kinetic energy release (Q value) of tritium decay is 
18.6 keV, only a very small fraction of decays falls into that region so that large statistics, very 
good background suppression, and an excellent energy resolution are required. If KATRIN should 
not see a positive signal, new experimental techniques would be required to push the sensitivity to 
even smaller m„. For example, it has been proposed to study nuclear recoils in bound state /3 decay 
of tritium 0] , to reconstruct the electron and nuclear kinematics in tritium decay j3] , or to measure 
the electron flux near the tritium endpoint in a storage ring 0]. However, all of these proposals 
are limited by the large Q value of tritium, which makes the neutrino mass a small effect. The 
decay 187 Re — > 187 0s offers a lower Q value of only 2.657 keV and thus better sensitivity to m u , 
but since it is a unique first forbidden decay, the small decay rate makes it difficult to accumulate 
sufficient statistics [9]. 
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Figure 1: An example for tuning of a Q value by ionization. Note that energy levels in the decay 194 Ir — > 
194 Pt* shown here have an uncertainty of 0(1.6 keV), so the figure is to be understood as an illustration of 
the principle idea only. 



In this paper, we investigate weak decays with even smaller Q values. In particular, we consider 
continuum /3 (c/3), bound state /3 (b/3), and electron capture (EC) decays. The key ideas are to 
consider decays to excited nuclear daughter states and to use ions instead of neutral atoms if 
necessary. As illustrated in Fig. [H an appropriate choice of the ionization level allows for some 
tuning of the Q value since every spectator electron contributes to Q with its energy gain or loss 



due to the change of the nuclear charge during the decay. For b/3 decay [lfj, lll|, ionization can 
also have the direct effect of opening up new decay modes. Our aim is to find decays that have 
sufficiently small Q values to depend appreciably on m u , but at the same time still have an absolute 
rate large enough to allow for a good signal-to-noise ratio. 

For b/3 and EC decay, the observable sensitive to m u is the decay rate, which could be measured 
by detecting gamma and x-ray photons accompanying the decay. In the case of c/3 decay, the 
sensitivity can be increased by guiding the decay electrons into a spectrometer similar to the ones 
used in Mainz, Troitsk, and KATRIN in order to also measure the /3 spectrum near the endpoint. 
For those decays where ionization is required to achieve sufficiently low Q, we propose to store the 
parent ions in a trap or in a storage ring. We discuss the feasibility of these ideas below. 



2. NUCLEAR DECAYS WITH ULTRA-LOW Q VALUES 

Nuclear decays with an ultra-small kinetic energy release Q <C 1 keV can occur only if the 
daughter nucleus has a state with excitation energy E* > fulfilling 

/3~: Qo - (B z +i,z - Bz,z) < E* < Q + Bz,i - -Bz+1,2 , 

P + : Qo ~ 2m e <E* <Q - 2m e + B z ,z - B Z -i,z , 

EC: Qo - B z ,i <E*<Q . 

Here, we have neglected contributions from the 0(eV) binding energies of outer shell electrons. Qo 
refers to the atomic mass difference of the parent and daughter nuclei 1 and Bz, n is the modulus of 
the total electron binding energy in an atom or ion with nuclear charge Z and n orbital electrons. 
In the first equation, (Bz+i,z — Bz,z) is the energy gain of the spectator electrons in the decay of 
a neutral atom. By ionization, the effective Q value can be reduced by up to that amount. In b/3 
decay, ionization can also increase Q by opening up decay modes to low-lying bound states. The 



1 Note that in most reference tables, the atomic mass difference is simply called Q, while we reserve that notation 
for the actual kinetic energy release in a decay, which is the quantity that determines the sensitivity to m v . 
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maximum possible increase occurs for b/3 decay of a hydrogen-like ion into a helium-like daughter 
state, and is consequently of 0(Bz,i — i?z+i,2)- Similarly, for j3 + decay, Q can be increased 
compared to (Qo — 2m e ), the value for neutral atoms, by removing spectator electrons and thus 
avoiding an energy loss of up to (Bz z — Bz—iz)- For EC, Q can be smaller than the atomic mass 
difference Qo by up to the binding energy of the Is electrons, which is of 0(Bz,\)- Q cannot be 
made significantly larger than Qo for EC. 

We list several candidate isotopes for low-Q cj3 ± , b/3~, and EC decay in Tabled! Since nuclear 
structure data is still very incomplete for many isotopes, it is quite possible that other suitable 
decays will be identified in the future. To keep the expected signal-to-background ratio large, we 
have only considered isotopes for which a low-Q decay is allowed from spin and parity arguments, 
while other decay modes (if present) are at least first forbidden or otherwise have a very small 
branching ratio. Note that, for some of our candidate isotopes, decay into the relevant excited 
daughter state E* has not been observed yet, so even though it is not forbidden by spin and parity 
arguments we cannot be sure that it exists. For each decay we have computed Q as a function of 
the electron configuration. The main uncertainties in this calculation come from the atomic mass 
differences Qo, which are typically known to O(keV) [12 . and from the binding energies of multi- 
electron configurations. We have estimated these binding energies using (I) the relativistic Hartree- 
Fock code atsp2K [H, 14] and (II) published atomic physics data and simulation results 151 IsT] 
(we only report the results of method (II)). The good agreement between the two independent 
estimates (I) and (II) shows that the atomic physics uncertainty in our Q values is < 100 eV and 
thus smaller than the uncertainties in most Qo values. An actual neutrino mass measurement 
would, however, require both, Qo an d the electron binding energies, to be known to an accuracy 
better than 0(m u ), and we discuss below how this could be achieved. Here, we deal with the 
uncertainties by reporting how small Q can be made if the present best fit values for Qo are taken 
at face value, and by how much Q can change if Qo is varied within present uncertainties. In all 
cases, we assume the degree of ionization and the daughter state E* to be chosen in the optimum 
way. 

For c/3 decay, the most promising isotopes at present are 188 W, 193 0s, and Ir with achievable 
Q values between and 1.3 keV, depending on the true value of Qo- A measure for the sensitivity 
of these low-Q c/3 decays to nonzero m u is the rate of events with electron energies in a small 
interval [Q — SE, Q] near the spectral endpoint. However, by considering the phase space factor 
and the Coulomb correction term (Fermi function) entering in the c/3 decay rate, it is easy to 
show that this number is independent of Q. To zeroth order, and neglecting differences in nuclear 
matrix elements, this seems to indicate that for achieving the same sensitivity as KATRIN in a 
low-Q experiment a similar number of stored parent atoms (10 19 ) would be required, which is far 



beyond the capabilities of present ion traps (< 10 6 -10 8 ) [19|,l2Gj] and storage rings (< 10 9 -10 11 ) [211 ] . 
However, the larger relative effect of m u makes a low-Q c/3 decay experiment more robust against 
many systematic errors. For instance, the required relative spectrometer resolution is smaller than 
in the tritium case. Moreover, if the time of each decay can be tagged by observing an associated 
gamma or x-ray photon, the spectrometer can be operated in the more sensitive time-of-flight 
(MAC-E TOFL mode [H]. Finally, it might be possible to combine our ideas with the methods 
proposed in [6-8] to measure also the energy and momentum of the recoil nucleus. All these effects 
should help to reduce the required number of stored ions, even though the experiment will still be 
extremely challenging. 

For b/3 decay, 163 Dy could provide Q ~ 1.5 keV. This isotope has the interesting property of 
being stable as a neutral atom, but becoming unstable to b/3 decay when ionized [22] • The most 
promising isotopes undergoing electron capture are 159 Dy and 163 Ho, for which M-capture with a 
very low Q value might occur even without ionization, depending on the exact value of Qo- 163 Ho 



has been studied previously in the context of calorimetric m u measurements in ref. [23J. In Fig. [21 
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Table I: Candidates for ultra-low Q decays. The calculation of Q values is based on data from 12, lBT- 18| 
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Figure 2: Relative effect of nonzero m v on the decay rate in b/3 and EC decay. 



we plot [T(m u = 0) — T(m u ^ 0)]/T(m u = 0)], the relative effect of nonzero m u on the b/3 or EC 
decay rate, as a function of Q. We see that even if Q ~ 100 eV is achieved, the effect of m u = 2 eV 
(0.2 eV) is only at the level of 10 -4 (10 -6 ). Even if all systematic uncertainties could be reduced 
to that level, detecting a deviation from the m u = case would still require the observation of 
few x 10 8 (10 12 ) low-Q decays. To complete the experiment within few years of measurement 
time, this would in turn require a very large and continuously replenished sample of about 10 16 
(10 20 ) stored parent particles 2 , implying extreme, and possibly prohibitive, requirements on isotope 
production and (in the case of ionized parent atoms) storage technology. Part of the problem is 
the fact that the nuclear matrix elements for the relevant decay mode are small. If they were of 



2 To arrive at this estimate, we have taken the known partial lifetime of the EC decay 159 Dy — > 159r Tb* with 
E* — 363.51 keV [l^] and have replaced the phase space factor and the electron wave function by the expressions 
appropriate for the low-Q decay. 
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0(1), the decay rate would be about 10 4 times larger. Let us emphasize again that decays with 
larger matrix elements (or even smaller Q values) may exist, but to identify them, more precise 
data on Qo values and on nuclear excitation levels is needed. 



3. FEASIBILITY OF A NEUTRINO MASS MEASUREMENT USING LOW-Q DECAYS 

In order to exploit the high m^-sensitivity of low-Q decays, one has to overcome several severe 
technological challenges. We consider the most important ones to be (A) producing a sufficient 
number of parent nuclei, (B) storing them, (C) obtaining an accurate prediction for the decay rate 
(for b(3 and EC decays, where no spectral information is available), and (D) counting the decays. 
In the following, we discuss some ideas on how these difficulties might be overcome. 

(A) Producing a sufficient number of parent nuclei. Most of the isotopes listed in 
Table U are unstable, so they would have to be produced artificially. At future facilities like FAIR 



at GSI, radioactive beams with at least 10 s - 10 10 ions/s can be produced [24| for nuclei not too 
far from stability. For isotopes with half lives of 0(days), this is in principle sufficient to sustain a 
sample of 10 13 -10 15 parent particles, but our discussion above shows that a competitive neutrino 
mass measurement would still require an improvement of several orders of magnitude unless a new, 
extremely favorable low-Q decay mode is discovered in the future. For an experiment using ionized 
parent atoms, an additional challenge is to remove ions in other than the desired charge state in 
order to avoid decays with larger Q value but identical experimental signature (i.e. identical 7 and 
x-ray fingerprint) as the considered low-Q decay. Due to the different charge-over-mass ratios of 
differently charged ions, this should in principle be possible. 

(B) Storing a sufficient number of parent particles. While an experiment using neutral 
atoms (e.g. 159 Dy and 163 Ho) can use a gaseous, liquid, or solid source, a setup using ions requires 
a trap or a storage ring. With present techn olog y, it is possible to store a total charge of 10 8 e 
(corresponding to 10 6 heavy ions) in a trap 1J^, 20] and 10 11 e in a storage ring 21]. Traps might be 



pushed to 10 9 e [25[ in the future, and the planned FAIR facility at GSI Darmstadt would provide 
storage rings with a capacity of 10 12 e [24(]. As already mentioned, this is still not sufficient to 
perform a low-Q (3 or EC decay experiment using ionized parent atoms competitive to KATRIN 
unless new decay modes with Q < 1 eV and large nuclear matrix element are discovered. This 
implies that, from the present perspective, decays of neutral atoms look more promising. 

(C) Predicting the decay rate T for b/3 or EC decay. The main unknowns in the 
computation of T are the nuclear matrix element, the nuclear mass difference, and the electron wave 
functions. To avoid the uncertainty in the matrix element, we propose to study not only the low-Q 
decay, but also a large-Q (i.e. high rate, but small m^-dependence) decay into the same nuclear 
final state to measure the nuclear matrix element. The mass difference Wq between the parent 
and daughter nuclei can be measured using ion trap mas s sp ectrometry. This technique currently 
provides an impressive relative accuracy of 0(1O -11 ) 26|, |27|, but for our purposes, this would still 



have to be increased by more than one order of magnitude to make the uncertainty in Wq smaller 
than the effect of the neutrino mass. The electron wave functions entering in T cannot be measured 
directly and have to be predicted by solving the multi-particle Dirac equation. The uncertainties of 
these predictions must be smaller than the expected effect of m u , but considering that many atomic 
x-ray spectra can be predicted to an accuracy below one per mille [2^], this could be feasible. To 
minimize the theoretical errors, one could 'calibrate' the numerical computation using experimental 
x-ray spectra, ionization energies, and other atomic physics data for the considered isotope. 

(D) Counting the number of decays. For b/3 and EC decay, the only observable sensitive 
to m u is the decay rate into the low-Q channel. To measure it, and to reject concurrent large- 
Q decay modes, we propose to detect characteristic gamma or x-ray photons accompanying the 
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decay. The main requirements for the photon detector are good solid angle coverage, high energy 
resolution, and efficient suppression of backgrounds from cosmic ray interaction products and 
radioactive impurities. To date, the best 7 detectors — employing extremely radiopure materials, 
active and passive shielding, and several meters of rock overburden — achieve background rates 
< 10 3 keV-V -1 and an energy resolution around 1 keV [29]. If the considered low-Q decay 
is accompanied by several photons, much better background suppression will be possible if the 
coincidence technique is used. Therefore, we estimate that backgrounds can be brought under 
control. 

4. CONCLUSIONS 

In this paper, we have discussed how continuum /?, bound state (3, and electron capture decays 
with extremely small Q values (<C 1 keV) can be realized and how they could be used to measure 
the absolute neutrino mass m u . To achieve sufficiently low Q values, i.e. sufficiently high sensitivity 
to m v , we have proposed to consider decays into excited nuclear daughter states, and, if necessary, 
to partially ionize the atoms to tune the electronic contribution to Q. We have discussed the tech- 
nological challenges that would have to be overcome in such an experiment, including production 
and storage of a large number of radioactive atoms or ions, obtaining accurate predictions for the 
decay rate as a function of m u , and counting the number of decays. We have found that the most 
promising decays to date are 159 Dy — > 159 Tb* and 163 Ho — > 163 Dy because, depending on the 
exact values of the respective atomic mass differences Qq, they may have low-Q EC decay modes 
even when neutral. Experiments using ions are much more challenging due to the large number of 
particles that must be stored. As a next step, it is crucial to measure precisely Qo for the isotopes 
listed in Table U in order to determine how small Q can be made for them. Also, more precise 
data on nuclear excitation spectra throughout the chart of nuclides is desirable in order to identify 
further candidates for low Q decays. 

Acknowledgments 

It is a pleasure to thank K. Blaum, F. Bosch, G. Heusser, B. Kayser, M. Lindner, Yu. Litvinov, 
Yu. Novikov, and M. Weber for useful and inspiring discussions. This work has been supported by 
the DFG-SFB TR 27 'Neutrinos and beyond'. Fermilab is operated by Fermi Research Alliance, 
LLC under Contract No. DE-AC02-07CH11359 with the US Department of Energy. 



[1] O. Elgaroy and O. Lahav, New J. Phys. 7, 61 (2005), hep-ph/0412075. 

[2] F. T. Avignone, III, S. R. Elliott, and J. Engel, Rev. Mod. Phys. 80, 481 (2008), 0708.1033. 

[3] C. Kraus et al., Eur. Phys. J. C40, 447 (2005), hep-ex/0412056. 

[4] A. I. Belesev et al., Phys. Lett. B350, 263 (1995). 

[5] J. Angrik et al. (KATRIN) (2004), available from http://bibliothek.fzk.de/zb/abstracts/7090.htm 

[6] S. G. Cohen, D. E. Murnick, and R. S. Raghavan, Hyperfme Interactions 33, 1 (1987). 

[7] M. Jerkins, J. R. Klein, J. H. Majors, and M. G. Raizen (2009), 0901.3111. 

[8] M. Lindroos, B. McElrath, C. Orme, and T. Schwetz, Eur. Phys. J. C64, 549 (2009), 0904.1089. 

[9] A. Monfardini et al., Prog. Part. Nucl. Phys. 57, 68 (2006), hep-ex/0509038. 

[10] J. N. Bahcall, Phys. Rev. 124, 495 (1961). 

[11] M. Jung et al., Phys. Rev. Lett. 69, 2164 (1992). 

[12] R. Firestone, Table of isotopes, http://ie.lbl.gov/toibook.html 



7 



[13] C. Froese Fischer, T. Brage, and P. Jonsson, Computational atomic structure: an MCHF approach 

(IOP Publishing Ltd., Bristol, UK, UK, 1997), ISBN 0-7503-0466-9. 
[14] C. Froese Fischer et al., The atsp2k code for atomic structure calculations, 

http://atoms.vuse.vanderbilt.edu/Elements/CompMeth/CompJVlethods.html. 
[15] G. Rodrigues, P. Indelicato, J. Santos, P. Patt, and F. Parente, Atomic Data and Nuclear Data Tables 

86, 117 (2004), ISSN 0092-640X. 
[16] W. R. Johnson and G. Soff, Atomic Data and Nuclear Data Tables 33, 405 (1985), ISSN 0092-640X. 
[17] D. R. Plante, W. R. Johnson, and J. Sapirstein, Phys. Rev. A 49, 3519 (1994). 
[18] J. A. Bearden and A. F. Burr, Rev. Mod. Phys. 39, 125 (1967). 

[19] F. Ames et al., Nucl. Instrum. Meth. A538, 17 (2005), ISSN 0168-9002. 

[20] I. Podadera Aliseda (2006), CERN-THESIS-2006-034, available from |http: / /cdsweb.cern.ch/record/ 
975263?ln=en. 

[21] N. Kalantar-Nayestanaki et al., Int. J. Mod. Phys. E18, 524 (2009). 

[22] K. Takahashi and K. Yokoi, Nuclear Physics A 404, 578 (1983), ISSN 0375-9474. 

[23] A. De Rujula and M. Lusignoli, Phys. Lett. B118, 429 (1 982). 

[24] GSI (2001), http://www.gsi.de/GSI-Future/cdr/, see also http://www.gsi.de/fair/index_e.html 
[25] K. Blaum, private communication. 
[26] K. Blaum, Phys. Rep. 425, 1 (2006). 

[27] K. Blaum, Y. N. Novikov, and G. Werth (2009), 0909.1095. 
[28] R. D. Deslattes et al., Rev. Mod. Phys. 75, 35 (2003). 
[29] G. Heusser, private communication. 



